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Abstract 

The Kac model describes the local evolution of a gas of N particles with three dimensional 
velocities by a random walk in which the steps correspond to binary collisions that conserve 
momentum as well as energy. The state space of this walk is a sphere of dimension 3A^ — 4. The 
Kac conjecture concerns the spectral gap in the one step transition operator Q for this walk. In 
this paper, we compute the exact spectral gap. 

As in previous work by Carlen, Carvalho and Loss where a lower bound on the spectral gap 
was proved, we use a method that relates the spectral properties of Q to the spectral properties 
of a simpler operator P, which is simply an average of certain non-commuting projections. The 
new feature is that we show how to use a knowledge of certain eigenfunctions and eigenvalues 
of P to determine spectral properties of Q, instead of simply using the spectral gap for P to 
bound the spectral gap for Q, inductively in N, as in previous work. The methods developed 
here can be applied to many other high-dimensional stochastic process, as we shall explain. 

We also use some deep results on Jacobi polynomials to obtain the required spectral informa- 
tion on P, and we show how the identity through which Jacobi polynomials enter our problem 
may be used to obtain new bounds on Jacobi polynomials. 
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1 The Markov transition operator Q for the Kac walk 

Let X]\f be the N particle state space consisting of A^-tuples v = {vi, . . . ,vn) of vectors Vj in 
with 

N N 

I P = 1 and Vj = . 

j=i j=i 
We think of a point v as specifying the velocities of N particles, and shall consider a random walk 
on Xtv that was introduced by Mark Kac [7j. At each step of this random walk, v is updated due 
to the effect of a binary collision that conserves energy and momentum — hence the constraints 
defining Xn. 

To specify this walk in more detail, we consider a collision in which particles i and j collide. 
Suppose that v* and v* are the post-collisional velocities, while Vi and Vj are the pre-collisional 
velocities. Then by momentum conservation, the center of mass velocity is conserved; i.e., 

V* + V* = Vi + Vj . 

Furthermore, by energy conservation, i.e., I^^*!^ + I'Vj P = + l^jP) ^-nd the parallelogram law, it 
follows that 

This leads to a natural parameterization of all the possible binary collision outcomes that conserve 
energy and momentum: The parameter o" is a unit vector in S"^, and when particles i and j collide, 
one updates v ^ v* = Rij^„{v) where 

Vi + Vj \Vi-Vj\ 

Vj = H —a 

2 2 

^ Vi + Vj \Vi-Vj\ 

Vi = —a 

^ 2 2 

vl = Vk for k^i,j . 

(1.1) 

The Kac walk on Xj\f is a random walk in which the steps are such binary collisions between 
pairs of particles. At each step, one picks a pair i < j uniformly at random, and also a unit 

vector a in S'^. One then makes the update described in (II. ip . Of course it remains to specify the 
probabilistic rule according to which a should be selected. In the physics being modeled here, the 
likelihood of selecting a particular a will depend only on the resulting scattering angle 9 , which is 
the angle between v* — Vj and Vi — Vj. In the parameterization above, this is the angle between a 
and Vi — Vj. That is, 

cos(fc') = a ■ 



Vi - V 



'3\ 



The scattering rate function 6 is a non negative integrable function on [—1, 1] with 



1 

b{u)du 

■1 
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Then for any Vi 7^ vj, and with da being the uniform probabihty measure on S^, 

b(a-^^^^)da = l . (1.2) 

(If = Vj, the colhsion has no effect, and can be ignored.) One selects a according to the probability 
density that is integrated in (jl.2p . 

There are several choices of b of particular interest. One is the uniform redirection model, given 
by b{x) = 1 for all —1 < x < 1. In this case, the new direction of the relative velocity, a is chosen 
uniformly from 5^. 

Another is the Morgenstern model |llj , [12] , or the uniform reflection model: For any unit vector 
uj ^ S"^, let be the reflection given by 

Hu){v) = V — 2{v • uj)uj . 

In the uniform reflection model, one updates the relative velocity according to 

Vi - Vj H^{vi - Vj) = V* - Vj 

with Lj chosen uniformly. The relation between to and a is given by a = H^{{yi — Vj) /\vi — Vj\), 
and computing the Jacobian of the map cj — > a, one finds 

h{x) 



^/2^/r^ ■ 

Both of these belong to the one parameter family 

6a(x) = (l-a)2"(l-x)^° . (1.3) 

Leaving the particular choice of h open, this completes the specification of the steps in the Kac 
walk. For more detail and background, see [7] and |3]. 

The main object of study here is the spectrum of the one step transition operator Q for this 
random walk, and the manner in which this spectrum depends on A^. Q is defined as follows: Let 
Vn be state of the process after the nth step. The one step Markov transition operator Q is given 
by taking the conditional expectation 

Q(t){v) = E((/)(iT„+i) \vn = v) , 

for any continuous function (j) on X^. 

From the above description, one deduces the formula 



N ^ ^ 



I m,jAv))b (<T ■ ^^^-^) da . (1.4) 
'^■Js^ \ \vi-Vj\J 



Let cjat denote the uniform probability measure on Xfq, which is the normalized measure induced 
on Xf^ as a manifold embedded in M^^. 

For any two unit vectors a and lo, one sees from (II. ip that 



CGL May, 2007 



4 



Prom this and the fact that the measure da at is invariant under v i-^ Rij^aV, it follows that for any 
continuous functions (j) and ip on Xj^, 



i{;{v)Q(j){v)daN = / / 'ip{Rij^ujv)(l){Rij,av)b{uj ■ a)dujdadaN ■ 
Jxn Js'^ Js^ 

It follows that Q is a self adjoint Markov operator on L'^{Xj\f, (Tat). Moreover, it is clearly a Markov 
operator; that is, in addition to being self adjoint, Q is positivity preserving and Ql = 1. 

The motivation for considering the spectral properties of Q stems from a theorem of Kac [7] 
that relates the continuous time version of the Kac walk to the nonlinear Boltzmann equation. 
For the details, see [7] or [3]. Let v{t) denote the random variable giving the state of the system 
at time t for the process run in continuous time with the jumps taking place in a Poisson stream 
with the mean time between jumps being Then the equation describing the evolution of the 

probability law of v{t), is called the Kac Master Equation: If the initial law on Xj\f has a density 
Fq, then the law at time t has a density F{v,t) satisfying 

^F{v, t) = N{I - Q)F{v, t) with F{v, 0) = Foiv) . 
ot 

The solution F{v, t) is of course given by 

F{v,t)=e''^Fo{v) , 

where 

£ = N{Q - I) . 

Since Q is a self adjoint Markov operator, its spectrum lies in the interval [—1,1], and since 
Ql = 1, the constant function is an eigenfunction with eigenvalue 1. It is easily seen that as long 
as b{x) is strictly positive on a neighborhood of x = 1, the eigenvalue 1 of Q has multiplicity one. 
It then follows that the spectrum of C lies in [— 2A^, 0], and that is an eigenvalue of multiplicity 
one. We impose this assumption on b throughout what follows. 

The Kac conjecture for this stochastic process pertains to the spectral gap 

a2 /„-r\ 



= inf \ — 4){v)C<f){v)d(JN 

Ixn 



(p {v)da]\f = 1 , / (j){v)da]\f = 
Xjv ■J Xff 



and states that 



lim inf Aat > . 



This was proved by Carlen, Carvalho and Loss [3j, but without an explicit lower bound. Kac also 
made a similar conjecture for a simplified model with one dimensional velocities and no conservation 
of momentum. Por this model, the conjecture was first proved by Janvresse [.6J, though her approach 
provided no explicit lower bound. The sharp bound for the simplified model was first established 
in [2]. See Maslen [9] for a representation theoretic approach. 

The main goal of the present paper is to compute exactly liminfjv^oo ^Af- We shall be able to 
do this under an easily checked condition relating A2 and the quantities 

If',., , „ I 



Bi = — xb{x)dx and B2 = — x b{x)dx . (1-5) 
2 J-i 2 J_i 

The condition, given in (jl.6p below, will turn out to be satisfied when b is given by ba, as in (jl.3|) . 
for all < a < 7/9. 
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1.1 THEOREM. Suppose that B2 > Bi and that 



A2 > —(1-^2) . 



(1.6) 



Then for all N >3, 




(1.7) 



Moreover, the eigenspace is three dimensional, and is spanned by the functions 



N 




a = 1,2,3 , 



(1.8) 



where v'j denotes the ath component ofvj. 



As we shall see in the next section, for many choices of b, including the ba, there is a simple 
monotonicity of the eigenvalues of Q for N = 2 which ensures that the eigenfunction providing the 
gap comes from a first degree polynomial, and thus that 



When ([Tg]) is satisfied, the condition ([L6|) reduces to (1 - -Bi)/(1 - S2) > 20/9. 

Next, notice that the eigenfunctions listed in (jl.Sp are symmetric under permutation of the 
particle indices. Indeed, the operator Q commutes with such permutations, so that the subspace 
of functions with this symmetry is invariant. As explained in [7] and [3j, it is the spectrum of Q 
on this subspace that is relevant for the study of the Boltzmann equation. 

Moreover, notice that in the collision rules (jl.ip . exchanging and v*^ has the same effect as 
changing a to —a. For this reason, if one's primary object of interest is the Boltzmann equation, 
one may freely assume that 6 is a symmetric function on [—1,1], since then replacing b{x) by 
(6(x) + b{—x))/2 will have no effect on the spectrum of Q on the invariant subspace of symmetric 
functions, or on the corresponding Boltzmann equation. (See the introduction of [4J for more 
discussion of this point in the context of the Boltzmann equation.) If B is symmetric, then Bi = 0, 
and we do have Bi> B2. 

However, it is interesting that the Kac conjecture holds without restriction to the symmetric 
subspace, and the that methods developed here can be used to determine the spectral gap even when 
b is not symmetric, and the eigenfunctions corresponding to the gap eigenvalue are not symmetric. 

When b is not symmetric, it may happen that Bi < B2. We shall give examples of this below. 
The next theorem gives the spectral gap and the eigenfunctions whenever A2 = 2(1 — Si), regardless 
of whether Bi < B2 or B2 < Bi. However, it gives the exact value of Ajv only for N > 7. Since 
we are interested in large values of N, this is fully satisfactory. Indeed, it is remarkable that the 
two theorems show that already at relatively small values of N, the behavior of the system is very 
close, qualitatively and quantitatively to the behavior in the large N limit. 



A2 = 2(l-Si) . 



(1.9) 
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1.2 THEOREM. Suppose that As = 2(1 - Bi). Then for all N>7, 

An = mm{ (1 - , (1 - B2) } • (1-10) 

Moreover, if B2 > Bi, the eigenspace is three dimensional, and is spanned by the functions 

N 

^(v) = ^\vj\\f a = 1,2,3, (1.11) 
i=i 

where v'j denotes the ath component ofvj. 

On the other hand, if B2 < Bi, the eigenspace is spanned by the functions of the form 

Ivil"^ - \vj\'^ and vf - vf , a = 1,2,3, (1.12) 

for all i < j. 

Finally, if Bi = B2, the eigenspace is spanned by both of the sets of functions listed in il.8]) 
and U.12\) together. 

For the family of collision rates introduced so far, the ba, one may apply Theorem 11.1^ as we 
have indicated, but only for a < 7/9. As we shall see in Section 2, Theorem 11.21 applies for all 
< a < 1, and in this case gives Ajy = {N/N — 1)(1 — B2) for N > 7. In order to illustrate the 
case in which Theorem 11.21 yields the gap Ajv = (N/N — 1)(1 — Bi), we introduce 

6„(x) = 2(a + l)l[o,i](x)x" a>0. (1.13) 

Since < x on (0, 1), it is clear that B2 < Bi for all a in this case. We show at the end of Section 
2 that at least for < a < 1, A2 = 2(1 — i?i), so that Theorem 11.21 applies in in these cases. 

The method of proof is quite robust, and in Section 10 we shall describe how it may be extended 
to determine the spectral gap of Q for still other choices of b that are not covered by the Theroems 
OandO 

The method of proof of these theorems relies on a basic strategy introduced in [3j, but which 
is extended significantly here. The strategy consists of exploiting an inductive link between the 
spectral gap of Q and the one of an operator P, an average over projections introduced in Section 
3. In fact, 

N 

An>j^—^{1-Iin)An-i (1.14) 

where 1 — is the gap of P. The eigenvalues of P are much easier to compute than the ones of 
Q since the range of P consists of sums of functions of single variables Vj. 

In the case of the original model treated by Kac, one is in the happy circumstance that Q has 
a single gap eigenfunction (f) which is also the gap eigenfunction of P for all A^, and when this is 
used as a trial function in the derivation of (jl.l4p . one sees that (jl.l4p actually holds with equality, 
giving an identity relating A^v and Aiv_i. Thus, starting at A = 2, where the gap can be easily 
calculated, the above formula yields a lower bound on Ajv that turns out to be exact. The model 
treated in this paper does not have this simplifying feature, even when the gap eigenfunctions of 
Q are also the gap eigenfunctions of P. Nevertheless, the ideas that lead to (jl.l4p can be used in 
such a way that we can still calculate the gap of Q exactly. Very briefly, here is how: 
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Let /i^ < fiN he any number and assume that there are finitely many eigenvalues /i^ < /i]^ < 
• • • < fJ'N^ < IJ-N P ■ Denote the corresponding eigenspaces by Ej. Let Vj be the smallest invariant 
subspace of Q that contains Ej. Lemma 4.1 in Section 4 provides us with the following dichotomy: 
Either 

N 

^N>^^{l-^^*N)^N-l (1.15) 

or else: 

The gap of Q is the same as the gap of Q restricted to ©jLi Vj. (1-16) 

If the threshold has been chosen so that the lower bound on Ajv provided by ()1.15p is at least as 
large as the upper bound on A at provided by some trial function in (B'jLiVj, then Aat is the gap of 
Q restricted to ©JLj^V^-. As we shall see, the Vj are finite dimensional, so determining the gap of Q 
on ®JL]^V^- is a tractable problem. In this case we have determined the exact value of A^r. 

To proceed to the determination of A^v for all large A^, one needs a strategy for choosing the 
threshold /i^. The lower the value of /u^ that is chosen, the stronger the bound (jl.lSp will be, but 
also the higher the value of m will be. The basis for the choice of is a trial function calculation, 
providing a guess A^r for the value of A^r. Indeed, natural trial functions can often be chosen on 
the basis of physical considerations. (The spectrum of the linearized Boltzmann equation is the 
source in the case at hand.) To show that the guess is correct, so that Atv = A at, we are led to 
choose jj,*^ so that 

~ N 

A^<^^(l-/i^,)A;v-i (1.17) 

Since A^v-i > A^v-i, this forces us into the second alternative in the dichotomy discussed above, 
so that the gap eigenfunction for particles lies in ®^^iVj. Indeed, if the physical intuition behind 
the guess was correct, the trial function leading to Aat will lie in ©^^V^-, and yield the gap. 

Choosing /i^ small enough that ()1.17p is satisfied might in principle lead to a value of m that 
depends on A^. However, in the case at hand, we are fortunate, and can work with a choice of 
that leads to a fixed and small value of m, but for which (I1.17P is satisfied for all sufficiently large 
values of A^ - hence the restriction to A^ > 7 in Theorem ll.2l 

As will be clear from this summary of the strategy, the determination of the spectrum of P is the 
main technical step that must be accomplished. As we mentioned before, this is relatively simple, 
compared to the determination of the spectrum of Q, since the range of P consists of functions 
that are a sum of functions of a single variable. 

For this reason, we can reduce the study of the spectrum of P to that of a much simpler Markov 
operator K acting on functions on the unit ball B in M^. In the analysis of K, we shall draw on 
some deep results on Jacobi polynomials [8], [13]. In fact, it turns out that the connection between 
our eigenvalue problems and pointwise bounds on Jacobi polynomials is through a simple identity, 
and applications of this identity can be made in both directions: We not only use bounds on Jacobi 
polynomials to bound eigenvalues, we shall use simple eigenvalue estimates to sharpen certain best 
known bounds on Jacobi polynomials, as we briefly discuss in Section 11. 

First however, we deal with a simpler technical problem, the computation of the spectral gap 
of Q for A^ = 2. 
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2 The spectral gap for N = 2 

For N = 2, the state space X2 consists of pairs {v, —v) with w G satisfying \v\'^ 
that for = 2 the colhsion rules (jl.ip reduce to 



1/2. Note 



a/V2 



and 



-a/V2 , 



since vi + V2 = 0- 

The map (f , — 1>) 1-^ V^v identifies X2 with the unit sphere S"^, and the measure da2 on X2 
with da on S^. Thus, we may think of Q as operating on functions on S^. In this representation, 
we have the formula 

Q(j){u) = I (f){a)b{u ■ a)da . 
Notice that if R is any rotation of M'^ 

{Q(j)){Ru)=l (t>{a)b{Ru-a)da= I (l){Ra)b{Ru ■ Ra)da = [ cl){Ra)b{u ■ a)d^ 
Js^ Js^ Js^ 



a 



oR){u) . 



That is, {Q4>) o R = Q[(j) o i?), and this means that for each n, the space of spherical harmonics of 
degree n is an invariant subspace of Q, contained in an eigenspace of Q. In turn, this means that 
we can determine the spectrum of Q by computing its action on the zonal spherical harmonics, i.e., 
those of the form Pn{e • u) where e is any fixed unit vector in M'^, and P„ is the nth degree Legendre 
polynomial. Now, for any function </>(u) of the form (t){u) = f{e ■ u), 



/ 4>{a ■ e)b{a ■ u)da . 



We choose coordinates in which u and e span the x, z plane with u 



' " 









and e = 





1 




t 



so that t = u ■ e. Then with a 



Qfit) 



sm U sm If 
cos 6 sin ip 
cos 9 



u) = Qf{e ■ u) where 



n f2n 



4vr./o JO 



f{t COS 9 + \/l — sin 9 cos ip)b{cos 0) sm9d0dLp 
-f I fits + \/l - tVl - cos ip)b(s)dsdip . 

tTTJo J-1 



(2.1) 



Now, if / is any eigenfunction of Q with Qf = A/, then evaluating both sides at t = 1, we have 



1 r 

A/(l) = - / /(s)6(s)ds. Thus, the eigenvalue is given by 



2i-i /(I) 
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As we have observed above, the eigenfunctions of Q are the Legendre polynomials. Thus, if P„ is 
the Legendre polynomial of nth degree with the standard normalization Pn(l) = Ij and A„ is the 
corresponding eigenvalue, 

1 

An = / Pn{s)b{s)ds . (2.2) 



2 7-1 

This explicit formula enables one to easily compute A2. For example, we can now easily prove 
the following: 

2.1 LEMMA. When b{x) = ba{x), as in ^.3\} . then 1 — B2 < I — Bi for all a < I, and moreover, 

A2 = 2(l-Ai) = %^ = (l-i?2)(3-a) , (2.3) 
z — a 

so that il.6\} is satisfied for all a with < a < 7/9. 
Proof: Using Rodrigues' formula 

1 d" 

pjx) = — -r-r-(a;^ - 1) 
and integration by parts, one computes 

{a)n {a)n 



A„ = (1 - a) 



(1 - a)n+i (2 - a)n 



where (a)„ = a{a + l)(a + 2) • • • (a + n — 1). Notice that for all < a < 1, decreases as n 
increases, so with the collision rate given by ha, 



A2 = 2(1 - Ai) = . (2.4) 

2 — a 



Next, one computes 



2(1 -a) , 4(1 -a) 

1-Bi = -i f and 1 - ^2 - 



a) (2 — a)(3 — a) 

Since 2 > 4/(3 — a) for a < 1, 1 — B2 < 1 — Bi for all a < 1. Moreover, from this computation, 
one readily obtains (|2.3|) and the statement concerning p.6p . □ 

In particular, the condition ()1.6p is satisfied in both the uniform redirection model (a = 0) and 
the Morgenstern model (a = 1/2). Thus in these cases we have the exact spectral gaps 

2 

A AT = for the uniform redirection model 

3 A - 1 



8 A 

A AT = for the Morgenstern model 

15 A - 1 



(2.5) 



We close this section with a remark that may provide a useful perspective on what follows. In 
determining the spectral gap of Q for N = 2, general symmetry conditions told us right away what 
all of the eigenfunctions were. A less obvious, though still simple, argument then provided us with 
the explicit formula (j2.2p for all of the eigenvalues. There is one last hurdle to cross: There are 
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infinitely many eigenvalues given by (j2.2p . and for a general b, we cannot determine which is the 
second largest by computing them all explicitly. What was particularly nice about is that in 
this case, the eigenvalues of Q were monotone decreasing: 

For other choices of b, this need not be the case. However, there are ways to use pointwise bounds 
on Legendre polynomials to reduce the problem of determining A2 to the computation of a finite 
number of eigenvalues using ()2.2p . For example, one has the classical bound (see Theorem 7.3.3 in 

my- 

\Pn{x)\' < —^=L= . (2.6) 

As long as b{x){l — is integrable, this gives an upper bound on A„ that is proportional to 

n~^/^: Define 

/ 1 \i/2 .1 
A„=(— J j b{x){l-x^)-'/^dx . 

Then, let no be the least value of n such that A^ < Ai. Then the second largest eigenvalue of Q is 

max An . 

l<n<rao 

We illustrate this by showing that for the rate function ba introduced in (??), A2 = 2(1 — Bi) 
at least for < a < 1. (Of course, the integrals in (j2.2p can be computed exactly in this case; see 
7.231, page 822 in [5]. however, we prefer to illustrate the use of (|2.6p ). 

By (IMI) and (HJ), 



/ /■! \ 1/2 / .1 X 1/2 

|An| < (a + l)f / x^^dxj { Pnixfdx 
a + 1 1 



< 



V2a + 1 ■ 



(2.7) 



Also, by ()2.2|) . Ai = i?i = (a + l)/(a + 2). Comparison of the formulas shows that for < a < 1, 
A„ < Ai for all n > 4. Thus it suffices to check that Aj < Ai for j = 2, 3 and 4 by direct computation 
with (j2.2p . Doing so, one finds that this is the case. Hence, Theorem 11.21 applies, and yields 
Atv = {N/N - 1)(1 - Bi) for > 7. 

Further calculation would extend this result to higher values of a. Notice that as a tends to 
infinity, ba{x) is more and more concentrated at x = 1, which corresponds to = 0. Thus, for large 
values of o, b(x represents a "grazing collision model". 

For N > 2, the operator Q is much more complicated, and direct determination of the spectrum 
is not feasible. Instead, we use an inductive procedure involving a auxiliary operator that we now 
introduce. 
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3 The average of projections operator P, and its relation to Q 

A simple convexity argument shows that for each j, 

sup{\vj\'^ : V £ Xn} = ^ ^ ^ ■ 

For each j, define TTj{v) by 

... rN' 



N-1 ' ' 

SO that TTj maps Xj\[ onto the unit ball B in M'^. 

For any function (j) in L^(Xjv, dcjv), and any j with 1 < j < A^, define Pj{(p) to be the orthogonal 
projection of (/) onto the subspace of L'^{X]\f,daN) consisting of square integrable functions that 
depend on v through Vj alone. That is, Pj{4>) is the unique element of L'^{Xj\f,da]\f) of the form 
f{-Kj{v)) such that 

/ 4>{v)g{'n-j{v))daN = f{TTj{v))g{lTj{v))d(7N 

for all continuous functions g on B. 

The average of projections operator P is then defined through 

p = ±yp. . 

If the individual projections Pj all commuted with one another, then the spectrum of P would 
be very simple: The eigenvalues of each Pj are and 1. Moreover, Pj4> = (p if ^^d only if (p depends 
only on Vj so that it cannot then also satisfy Pk4> = (j) foi' ^ 7^ Ji unless (p is constant. It would then 
follow that the eigenvalues of P would be 0, 1/A^ and 1, with the last having multiplicity one. 

However, the individual projections Pj do not commute with one another, due to the nature of 
the constraints defining Xjy. 

We now define 



;UAr = sup<^ / 4>{v)P(t){v)daN / cl){v)daN = l, / 4>{v)daN = > ■ (3.1) 

The P operator is simpler than the Q operator in that if (p is any eigenfunction of P with non-zero 
eigenvalue, then clearly (p has the form 

N 

i=i 

for some functions fi, ■ ■ ■ , Jn on B. For > 4, most of the eigenfunctions of Q have a more 
complicated structure. Nonetheless, there is a close relation between the spectra of Q and P, as 
we now explain. 

To do this, we need a more explicit formula for P, such as the formula (jl.4p that we have for 
Q. The key to computing Pj(p is a factorization formula [3] for the measure a^- Define a map 
T]\r : X]y^i X B ^ X^ as follows: 



TN{y^v)=[a{v)yi--j^=v,...,a{v)yN-i-^^^=^v,^^^v\ , (3.2) 
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where 



a'^(v) = 1 - luP . 



This map induces coordinates {y,v) on X]\f, and in terms of these coordinates, one has the 
integral factorization formula 



1 53^- 



|^3Af-4| 



B 



(j)(Tiy{y,v))daN- 



(1 _ |^;|2)(3^-8)/2d^ . 



It follows from this and the definition of P/v that 

Pn4>{v) = f o TTNiv) 

where 



fiv) 



(j)(TN{y,v))daN- 



For j < N, one has analogous formulas for Tj and Pj, except the roles of vn and Vj are 
interchanged. 

Next, we make the definition for Q that is analogous to (j3.ip for P: Define A at by 



\n = sup <j / (l){v)Q(l){v)daN 
'Xn 



4> {v)da]s[ = 1 , / <l){v)da]\[ = 



(3.3) 



With this explicit formula in hand, and the definitions oi fiN and A at, we come to the funda- 
mental fact relating P and Q: 

3.1 LEMMA. For any square integrable function (j) on X]\f that is orthogonal to the constants, 

{^,Qct>) < Xn-iUWI + (1 - Xn-i){^,P(P) , (3.4) 
where (•,•) denotes the inner product on LP' {X i\i , a n) ■ 

Proof: To bound {(j), Qcp) in terms of Aat-i, define for 1 < k < N , the operator Q^^^ by 



N - 1 
2 



-1 



V / 0(i?„-,(^7))da . 



That is, we leave out collisions involving the fcth particle, and average over the rest. Clearly, 

fc=i 

Therefore, for any 4> in L?'{Xn ^gn), 

N 



k=i 
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Using the coordinates {y,v) induced by the map : Xjv-i x B ^ X]\f, it is easy to see that for 
^ 7^ ^) J 7^ ^1 Ri,j,<j ^cts only on the y variable. That is, for such i and j, 

Ri,jA'^k{y,v)) =Tk{Ri,j,a{y),v) . 

Thus, if we hold v fixed as a parameter, we can think of {Q^''^(/)){Tk{y, v)) as resulting from applying 
the — 1 dimensional version oi Q to (p with held fixed. 

To estimate Atv, we need estimate {^,Q<j)) when (j) is orthogonal to the constants. When (p is 
orthogonal to the constants, and we fix v, the function 

y ^ <P{Tk{y,v)) 

is not, in general, orthogonal to the constants on Xtv-i- However, we can correct for that by adding 
and subtracting Pk4>- Therefore 

{{(/)- Pk<P),Q^''\<P-Pk(l))) < XN-iU-Pk^l 

= ^NM\mi-\\Pk(t>\\i) 



Ajv-i(||0||i-(0,Pfc0)) . 



(3.5) 



Then since Q^^)Pk(t) = Pk(t) and since is orthogonal to — 

= m - Pk<i)) + Pk<i))Q^^\{(t) - Pk(i)) + Pk<i))) 

= {{^-Pk(l)),Q^''\(l)-Pk<P) + {Pk(l>,Pk<l>) 
= ((</.- P,.0),Q('=)(0-Pfc</.)) + (0,Pfc0) 

< XN-i{\ml - {4>Pk4>)) + {(t>Pk^) 

(3.6) 

Averaging over k, we have p.4p . □ 

Lemma l3.ll was used as follows in [3]: Any trial function (j) for Aat is a valid trial function for 
fJ-N, so that 

Aat < AAf_i + (1 - XN^i)fJ-N ■ (3.7) 
Then since A^v = N{1 — Aat), we have 

N 

An > - m)^N-i . (3.8) 

N 

Therefore, with un = — — -{^ — /^Af), for all > 3, 




Thus, one route to proving a lower bound on Aat is to prove an upper bound on fijy, and hence an 
lower bound on aj^. This route led to a sharp lower bound for A^v — the exact value — for the 
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one dimension Kac model investigated in p|. However, it would not lead to a proof of Theorem 
II. li The reasons for this are worth pointing out before we proceed: 

As we shall see below, the eigenspace of P with the eigenvalue fiN — the gap eigenspace of P 
— is spanned by the functions specified in (jl.Sp . Granted this, and granted Theorem ll.il whenever 
condition (jl.6p is satisfied: 

For (1 - B2) < (1 - Bi) , = Xn^ P(p = , 

while 

For (1 - Bi) < (1 - B2) , Q(l) = Xn(I) P(t>^ i^N(t> ■ 



In the second case, (1 — Bi) < (1 — B2), the mismatch between the gap eigenspaces for Q and 
P means that equality cannot hold in (j3.7p . and hence the recursive relation (j3.8p cannot possibly 
yield exact results in this case. 

Moreover, even in the first case, (1 — B2) < (1 — where there is a match between the gap 
eigenspaces of Q and P, there still will not be equality in (j3.7p . The reasons for this are more 
subtle: The inequality (|3.7p comes from the key estimate (|3.6j) . Considering (13. 6p . one sees that 
equality will hold there if and only if 

for each k, where (0 — Pk4') is regarded as a function on Xjv-i through the change of variables 
Tfc : (X7V-1, B) Xn that was introduced just before Lemma [3Tl 

However, if (j) is in the gap eigenspace for Q on X^, Theorem 11.11 tells us that it is a linear 
combination of the three functions specified in (jl.Sp . all of which are homogeneous of degree 3 in 
Because of the translation in (j3.2p . which is due to momentum conservation, [(f) — Pk4>) is regarded 
as a function on X7V-1 will not be homogeneous of degree 3 — it will contain lower order terms. 
Hence {(j) — Pk4>) will not be in the gap eigenspace for Q(^'). 

The main result of the next section provides a way to use more detailed spectral information 
about P to sharpen the recursive estimate so that we do obtain the exact results announced in 
Theorem 11.11 

4 How to use more detailed spectral information on P to deter- 
mine the spectral gap of Q 

The following lemma is the key to using (13. 4p to obtain sharp results for the model considered here. 
4.1 LEMMA. For any N >3, let fi^ be a number with 

such that there are only finitely eigenvalues of P between and fiN-' 

(tti) (1) 
< H}^ '<■■■< fl}^' < ■ 
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Let /i]Y denote and then for j = 0, . . . ,m, let Ej denote the eigenspace of P corresponding 
(i) 

to /x]^ . Let Vj denote the smallest invariant subspace of Q that contains Ej. Let Vj he the largest 
eigenvalue of Q on Vj. 
Then either 

Aat = max{fo, • • • , i^m} , (4.1) 

or else ^ 

^n>j:^{1-Ii*n)^n-i (4.2) 

If lJi\[ = IJ'^i^\ then we have the same alternative except with strict inequality in 

Proof: If Atv > max{z/o, . . . , t'mji then in the variational principle for Aat, we need only consider 
functions (p that are orthogonal to the constants, and also in the orthogonal complement of each of 
the Vj. This means also that (p belongs to the orthogonal complement of each of the Ej. But then 

Using this estimate in (j3.4p . we have (j4.2p . Moreover, if //^ = then strict inequality must 

hold in the last inequality. □ 

Lemma 14.11 gives us the dichotomy between (jl.lSp and (jl.l6p that plays a key role in the 
strategy described in the introduction. To put this strategy into effect, we must first carry out a 
more detailed investigation of the spectrum of P. The main result of the next section reduces the 
investigation of the spectrum of P to the study of simpler operator — the correlation operator K, 
which is a Markov operator on functions on the unit ball B in R^. 



5 The correlation operator K, and its relation to P 

While Q and P are both operators on spaces of functions of a large number of variables, the 
problem of computing the eigenvalues of P reduces to the problem of computing the eigenvalues of 
an operator on functions on B, the unit ball in M'^: 

First, define the measure i^n on B to be the "push forward" of aN under the map iTj. That is, 
for any continuous function f on B, 

fiv)di^N = / f{T^jiv))daN ■ 

B JXn 

By the permutation invariance of cjat, this definition does not depend on the choice of j. By direct 
calculation [3], one finds that 

d-A^(-) = P^(l-H')^'''''^/'d.. (5.1) 

Now define the self adjoint operator operator K on L^{B, dviy) through the following quadratic 
form: 



{f,Kf)L2^,)= /(7ri(t;))/(7r2(iT))daiv (5.2) 
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for all / in L'^{B,di'j^). Equivalently, 

(ir/)o7ri =Pi(/o7r2) . (5.3) 

Note that by the permutation invariance of crjv, one can replace the pair (1, 2) of indices by any other 
pair of distinct indices without affecting the operator K defined by (15. 3|) . This is the correlation 
operator. 

To see the relation between the spectra of P and the spectra of K, suppose that (f) is an 
eigenfunction of P that is symmetric under permutation of the particle indices. (These symmetric 
eigenfunctions are the ones that are significant in the physical application.) Then since any vector 
in the image of P has the form X^jLi fj ° for functions /i, . . . , /at on B, we must have, for 
symmetric, 

N 

i=i 

Now we ask: For which choices of / will (p given by ()5.4p be an eigenfunction of P? To answer 
this, note that by by (jS.Sp . 

N 

Pk(t> = f°'^k+ Yl PkU°^3) ■ (5.5) 

Therefore, from (|5.5p and the definition of K, P^cf) = /oTr^ + (A^ — l)(i<'/) ovr^. Thus, averaging 
over fc, 

N 

1 A'' — 1 

P<t>=J^<t^ + ^^Y.^Kf)o^^. (5.6) 

i=i 

In the case Kf = nf, this reduces to 

P(t> = ^{l + {N-1)K)4> , 

and thus eigenfunctions of K yield eigenfunctions of P. It turns out that all symmetric eigenfunc- 
tions arise in exactly this way, and that all eigenfunctions, symmetric or not, arise in a similar way, 
specified in the next lemma. 



5.1 LEMMA. Let V be the orthogonal complement in L'^{Xn,(Tn) of the kernel of P. There is a 
complete orthonormal basis of V consisting of eigenfunctions (p of P of one of the two forms: 

N 

(i) For some eigenfunction f of K, cp = f ° "^k- ihis case, if Kf = nf , then Pep = fi(p where 

k=l 

^ = ^ (1 + (iV - 1)k) . (5.7) 

(a) For some eigenfunction f of K , and some pair of indices i<j,(p = foTTi — fonj. In this 
case, if Kf = nf , then Pep = ^(p where 
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Proof: Suppose that is an eigenfunction of P with non zero eigenvalue fx, and (p is orthogonal 
to the constants. By the permutation invariance we may assume that either (p is invariant under 
permutations, or that there is some pair permutation, which we may as well take to be cri_2, such 
that (j) ° en, 2 = —0- We will treat these two cases separately. 

First suppose that <j) is symmetric. We have already observed that in this case, the recipe 
(j) = X^jLi / ° ^j; with / an eigenfunction of K, yields symmetric eigenfunctions of P. We now 
show that all symmetric eigenfunctions of P on y have this form. 

First, simply because such a </> is in the image of P, and is symmetric, seen that must have 
the form (j5.4|) . It remains to show that / must be an eigenfunction of K. Then by (j5.6p . ^i(f) = P(j) 
becomes 

N N 
llY.f0TTu = -Y.^f + {N-l)Kf)07,k . 
k=l k=l 

Apply Pi to both sides to obtain 

1 ([/ + {N- l)Kf] + {N- l)K [f + [N- l)Kf\) = + {N- l)Kf) 

which is 

1 (/ + (iV - l)Kf f = ^{I + {N- l)K)f . (5.9) 

Since /x 7^ 0, / is not in the null space of either I + {N - l)K or {I + {N - l)Kf . It then follows 
from ([521) that 

l(/ + (iV-l)if)/ = /./. 

Thus, when (j) is symmetric, there is an eigenfunction f oi K with eigenvalue k, such that (p = 
EfcLi / ° ^fe and 

^ = ^(l + (iV-l)^) . 
We next consider the case in which 

(j) o (Tl,2 = —4> ■ 

Note that 

Pk{<poai^2) = Pk(p = ^ 
whenever k is different from both 1 and 2. It follows that 

N 



k=l 



The right hand side is of the form f{vi) — f(v2), and hence (p must have this form if it is an 
eigenvector. Taking (p = f o m — f o tt2 we have 



1 ^ 1 

^Y.Pk<P = ^i{f-Kf)o7r,-U- Kf) o TTa) 



k=l 

Hence when Pep = fi(p and (p is antisymmetric as above, There is an eigenvalue k of K such that 

1 - K 
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This proves the second part. □ 

Lemma 15.11 reduces the computation of the spectrum of P to the computaton of the spectrum 
of K. We undertake this in the next three sections. 



6 Explicit form of the correlation operator K 

For any two functions / and g on B that are square integrable with respect to vn, consider the 
bihnear form /(vri(iT))(7(7r2(w))d(TAr. It is easily seen from ()5.3p that 



if, Kg) = / fi7riiv))g{Tr2iv))daN , 



where here, (•, •) is the inner product on L^{B, vn)- 

Computing the right hand side using the factorization formula (13. 2p . but for Ti instead of T/v, 
one finds, for N > 3: 

The explicit form of K is slightly different for = 3. We can see this different form as a limiting 
case, if we make the dimension a continuous fact. The following way of doing this will be convenient 
later on: 

For Q > — 1, define the constant Cq by 



C„= (^^(l-lyprdy) ' 



so that for 



3iV-8 

a = 



2 

2\«, 



and then 



N-1 ^ ' ' " A^-1 

Now, as approaches 3, a — 3/2 approaches —1. Then the measure Cq(1 — |yp)"dy concentrates 
more and more on the boundary of the ball B, so that in the limit, it becomes the uniform measure 
on 5^. Understood in this way, the formula remains valid at a = 1/2; i.e., at = 3. 

It is clear that is a self adjoint Markov operator on L'^{B,i'n), and that 1 is an eigenvalue 
of multiplicity one. With more effort, there is much more that can be said; the spectrum of K can 
be completely determined. 
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7 The spectrum of K and ratios of Jacobi polynomials 

In studying the spectrum of the correlation operator, it is in fact natural and useful to study a 
wider family of operators of this type. Fix any a > 1/2, and any numbers a and b such that 

+ 6^ = 1 . 

Then define the generalized correlation operator, still simply denoted by K, through 

Kgiv) = C,_3/2 9 (ayr^N^y + bv) (1 - \y\^r-'/^dy . (7.1) 

Notice that as v and y range over B, the maximum of loi^l — \v\'^y + bv\ occurs when ay and 
bv are parallel. In that case, 

\a^l-\v\^y + bv\ = |a||2/|Vl- |i;|2 + |6||^;|l < {a^ + b^)y\{l - \v\^)\y\^ + \v\^)^/^ < 1 . 

Thus, as V and y range over B, so does 

u{y, v) = a^/l — \v\'^y + bv , (7.2) 



and g{a^J 1 — \v\'^y + bv) is well defined for any function g on B. Thus, K is well defined. 
Now when 

"= N-1 ^"""^ ^ = ~iV^' ^^-^^ 

wc know that K is self adjoint because in that case it is defined in terms of a manifestly symmetric 
bilinear form. We shall show here that K is always self adjoint for all + b'^ = 1, and that the 
eigenvalues of K are given by an explicit formula involving ratios of Jacobi polynomials. 

To explain this, we fix some terminology and notation. For any numbers a > — 1 and /? > — 1, 
pjl°''^^ denotes the nth degree polynomial in the sequence of orthogonal polynomials on [—1, 1] for 
the measure 

{l-x)"{l + x)f^dx , 

and is referred to as the nth degree Jacobi polynomial for {a, (3). As is well known, {Pn'^'^^}n>o is 
a complete orthogonal basis for L'^{[—1, 1], (1 — a;)"(l + x)^dx). 

Of course, what we have said so far specifies pjf''^^ only up to a multiplicative constant. One 
common normalization is given by Rodrigues' formula 

(_-\\n r]'"' / \ 



For this normalization, 

p(-/^)(l)=f ) and P(-^)(-l)=( I . (7.4) 
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7.1 LEMMA. Fix any a > 1/2, and any numbers a and b such that + b"^ = 1, and define 
K through the formula (7.1). Then K is a self adjoint Markov operator, and the spectrum of K 
consists of eigenvalues Kn,e enumerated by non negative integers n and i, and these eigenvalues are 
given by the explicit formula 

where (3 = 1 + 1/2, a is the parameter a entering into the definition of K . 

Proof: To see that K is self adjoint, we write it as a bilinear form, and change variables to reveal 
the symmetry. The change of variable that we make is naturally (y, v) (n, v) with u{y, v) given 
by (|7.2p . From (|7.2p . one computes y = u — bv/{ay^l — |f P, so that 



i-\y\' 



— a^\v\^ — \u\^ — — 2bu ■ v 

62(1- |t,|2) 

b'^ - (\u\^ + \v\'^) - 2bu ■ V 



0^(1 — |t)P 

The Jacobian is easy to work out, and one finds that dudv = a^(l — \v\'^)^/'^dydv, so that 
f{v)Kg{v)Ca{l-\v\''rdv = 

f{v)g{u)a-^^ [a^ - {\u\^ + \v\^) - 2bu ■ v] C^^^/^dudv . 

B JB 



(7.6) 




(7.7) 



This shows that the operator K is self adjoint on L'^{B,Ca{^ — I'^P)") foi" all a > 1/2, and all a 
and b with 0^ + 6^ = 1. 

Our next goal is to prove the eigenvalue formula (17. 5p . This shall follow from several simple 
properties of K. 

First, K commutes with rotations in M'^. That is, if i? is a rotation on M^, it is evident that 

K{g o R) = (Kg) o R . 
Hence we may restrict our search for eigenfunctions g of to functions of the form 

g{v) = h{\v\)\v\'ye,rn(.v/\v\) 

for some function h on [0, 00), and some spherical harmonic 3^^,m- 

Second, for each n > 0, K preserves the space of polynomials of degree n. To see this notice 
that any monomial in y^l — \v\'^y that is of odd degree is annihilated when integrated against 

any even monomial in y^l — \v\'^y is a polynomial in v. 

Combining these two observations, we see that K has a complete basis of eigenfunctions of the 
form 

gn,e,m{v) = hn,i{\v\'^)\vfy£,miv/\v\) 

where /i„ ^ is a polynomial of degree n. 
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To determine these polynomials, we use the fact that K is self adjoint, so that the eigenfunctions 
9n,e,m can be taken to be orthogonal. In particular, for any two distinct positive integers n and p, 
the eigenfunctions gn,e,m and gp/^m are orthogonal in L'^(B,Ca{i — Hence for each £, and 

for n ^ p, 

/in,Kkl')/ip,Kkl')(l - IHTlH^'dT; = . 



'M<i 

Taking r = |vp as a new variable, we have 







This is the orthogonality relation for a family of Jacobi polynomials in one standard form, and this 
identifies the polynomials hn/- A more common standard form, and one that is used in the sources 
to which we shall refer, is obtained by the change of variable t = 2r — 1, so that the variable t ranges 
over the interval [—1, 1]. Then for a, f3 > —1, Pn"'^\t) is the nth degree orthogonal polynomial for 
the weight (1 — t)"(l + 1)^ . With the variables t and |f p related as above; i.e., 

t = 2\v\'^ - 1 , 

KA\v\') = Pt^\t) 

for 

Now that we have all of the eigenfunctions determined, a further observation gives us a simple 
formula for the eigenvalues. Consider any eigenfunction g with eigenvalue k, so that Kg{v) = Kg{v). 
Let e be any unit vector in R^. Then since g is a polynomial and hence continuous. 



lim Kg{te) = lim / g (a^/l^y + Me) 0^.3/2(1 " l2/|)"-'/'d 
= 9 {be) , 



y 



(7.8) 



since Kl = 1. Combining this with Kg{v) = K.g{v), we have 

g (be) = Kg{e) . 

Now consider any eigenfunction gn,i,m of the form given above, and let Hn,e be the corresponding 
eigenvalue, which will not depend on m. Then taking any e so that 3^£,m(e) 7^ 0, we have that 

'^n,e = T — 77rb ■ (7.9) 

Changing variables as above to express this as a ratio of Jacobi polynomials, we finally have proved 
(173]). □ 
One might expect the largest eigenvalues of K to correspond to eigenfunctions that are poly- 
nomials of low degree. After all, in a system of orthogonal polynomials, those with high degree 
will have many changes of sign, and one might expect considerable cancelation when applying an 
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averaging operator, such as K, to them. Therefore, let us compute the Kn,£ for low values of n and 
£. We find from ()7.5p . using the value b = —1/{N — 1) from (j7.3p . that 

Ko,i = Ki,o = ' (7-10) 

so that is negative for n + i = 1. For n + i = 2, we find from ()7.5p . 



5iV-3 



K2,o(A^) 



3(A^- 1)3 
(iV-3)(15A^2 _i5^_^4) 



Evidently, for large A^, 



3(3iV-4)(iV- 1)4 



(7.11) 



while 



Thus, one might expect that at least for large values of N, 1, ki^i(A^), K2fi{N) and ko,2(-^) are the 
four largest eigenvalues of K, and that kq,! = i^ifi is the most negative, with all other eigenvalues 
of K lying strictly between these. 

We shall show in the next section that this is indeed the case for all > 4, and that 1 and ki^i 
are the two largest eigenvalues of K for all > 3. 

When we use Lemma l5.ll to convert this to spectral information on P, we find that ko,1) ^i,o 
and Kg, 2 all correspond to the same eigenvalues of P, namely 

1 / 1 \ 1 / / X 1 \ 1 



N\ N-lJ A^V ^(A^-l)V A^-1' 

This is the eigenvalue of P that shall play the role of ^j^T^ iii our application of Lemma 14.11 

Let us conclude this section by recording a number of useful calculations that can be made 

using (j7.5p . 

For A^ = 3, we have 

1 13 1 

'^i.ilS) = 2 > ^2,2(3) = — > ^0,2(3) = - > K2,o(3) = . (7.12) 



For A^ = 4, we have 



For A^ = 5, we have 



17 1 23 

'^1,1(4) = - > /^o,2(4) = - > Ai2,o(4) = — . (7.13) 



Ki,i(5) = ^ > ^2,o(5) = ^ > '^0,2(5) = ^ . (7.14) 
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In each case, the second largest eigenvalue after 1, among the ones listed, is ki^i. In the next 
section we shall see that the list is not misleading: ki^i is the gap eigenvalue. However, note that 
the third largest eigenvalue comes from different values of n and £ for each of = 3, = 4 and 
N = 5. As we shall see, things do settle down for N > 5; the third largest eigenvalue does turn out 
to be K2,o in a-ll such cases. 

7.2 LEMMA. For all N >5, ki,i(A^) > K2,o{N) > ko,2(A^). 
Proof: From (ffTTIl . 

,.r^ 2N(3N^ -15N + 8) 



3(3Af-4)(A^- 1)4 ■ 

A simple calculation shows that the roots of the polynomial in the numerator are less than 5, so 
that H2,o{^) > '*o,2(A'^) for A^ > 5. A similar argument applied to ki.i(A^) — K2,o(-^) yields the 
conclusion of the lemma. □ 

Our goal in the next section is to show that for all A^ > 4, there are no eigenvalues Kn/ with 
n+i > 2 that are larger that the ones listed above, and that for A^ = 3, the three largest eigenvalues 
are 1 = ft;o,o > 1/2 = ki^i > 13/40 = ^2,2- However, since there is no simple monotonicity in n + ^, 
this shall require some detailed estimate on ratios of Jacobi polynomials. 

We shall also need to know that in all case ko,i = ki,o = ~l/(-^ ~ 1) is the most negative 
eigenvalue. This will tell us the four largest eigenvalues of P for A^ > 4, and the three largest for 
n = 3, and this shall turn out to be enough to prove the main result, Theorem ll.il 

Finally, the value of K2,2(A^) will play an important role in the proof of Theorem [L2l and so we 
record the expression here: 

,,,, 21Ar3 _ 50^2 ^ 27A^ - 4 ,^ ^ , 

"^^^l'"' = (3Ar-4)(iV-l)« • <"^' 



8 The determination of the spectrum of K 

The main result in this section is the following theorem: 

8.1 THEOREM. For N > 5 and all n and I with n + £>2, 

- < K^A^) < ^0,2 (A^) • 

For N = 4 and all n and i with n + £ > 2 

1 



< K„,£(4) < K2,o(4) . 



N -1 

For N = 3 and all n and £ with n + £ > 0, except for n = 1,£ = 1, 



1 13 

< UnA^) < ^2,2(3) = — . (8.3) 



iV - 1 - ' - "'"^ ' 40 
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We present the proof at the end of this section after a number of preparatory lemmas. These 
lemmas rest on two deep results about Jacobi polynomials. One is a formula due to Koornwinder 
[8] (see also |lj, pp. 31 jff.)that was already applied in [3]: 

For all — 1 < < 1, all n and all a > /3, 



1 + X — (1 — x)r 



+ i\/ 1 — x^r cos{9) 



dmaAr,0) (8.4) 



jr^il) Jo Jo 
where 

rriaA^, 9) = c^A^ - r2)"-/3-V2/3+i (sin^)^^ drdO , 

and Ca,i3 is a normalizing constant that makes dma,i3 a probability measure. 

Koornwinder's bound is very useful for obtaining uniform control in n for given i and N. But 
since in Lemma |7.H 

OAT C 1 

a = and ^ = £ + -, (8.5) 

we can only apply (j8.4p when 

^<r = '-^ (8.6) 

As in [3], one may use this formula to show: 

8.2 LEMMA. For all i with 2<l<t, and all n > and N >3, 

KAN)\ < (]v3l)2 = ''0'2(iV) . 

Note that while this lemma does not address the case n = 0, this is not a problem: we have the 
explicit formula 

«°''=(¥^)'- ('*^^> 

To handle large values of i, we need another deep result, which is a uniform pointwise bound 
on the orthonormal Jacobi polynomials that was obtained by Nevai, Erdelyi, and Magnus. [13j: 
Let pn'^ be the orthonormal Jacobi polynomial of degree n with positive leading coefficient for the 
weight w{x) = (1 — x)"(l + x)^. It was shown in [I3] that for all a > —1/2 and /? > —1/2 and all 
non negative integers n, 



2e(2 + ^/^TP^) 



^xG[-l,l] - a;2y;(x)p°'^(x)2 < 



vr 

Of course, we could use the orthonormal Jacobi polynomials in the ratio formula (j7.5p . since 
any normalization factor would cancel out in the ratio. However, the exact formula (j7.4p for the 
denominator in (j7.5p is simplest in the other normalization. Hence we need the relation between 
Pn'^ and Pn'^ , which is given by pn'^ = InPn'^ where 



1/2 



2n + a + P+lT{n + l)T{n + a + /? + 1) 
2^+^+1 T{n + a + l)T{n + (3 + 1) 

Therefore 

Pn'^(x)^ ^ 1 2er(n + ifTja + 1)^(2 + ^a^ + P^) 
P"'^(l)2 " II Vl^wix)TTr{n + a + 1)2 
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At this point it is perhaps worth noting that since the spectrum of K lies in [—1, 1], any upper 
bound on its eigenvalues by a number larger than one is vacuous. This implies that for certain 
regions the identity (j7.9p will provide a stronger bound than (8.9). We shall return to this point at 
the end of the paper. 

2 13 
Substituting x = —1 + ^-y^-, 1^ = ^ + —, a = — — 4 in (jS.Sp . and then multiplying by 

1 

yields 



{N - 1)2' 
where 

where 



~^IA^) = -9i{n,i,N)g2{N)g3{n,N)g^{n,£,N) (8.11) 



vr 



9iin,i,N) 

92iN) 
93{n,N) 
g4{n + l,N) 



U + \/9Af2 - 48iV + 65 + 4^2 + 
I 3A^ + 4n + 2^ - 5 



\N{N - 2) 
r(n + l)r(|A^-3) 

r (n + fiV - 3) 
{N - l)2r (n + £ + 1) r (|iV - 3) 

r{n + e + lN-l) 



^.12) 



Our goal now is to extract a reasonably tight upper bound for K,n/{N) with as as much mono- 
tonicity in n, i and as possible. The next lemmas address this goal. 

8.3 LEMMA. For £>0, N > 3, and n > 



gi(n,£,N)<{ + 1 , (8.13) 

' ' ^ - * 3A^ + 4n + 2£-5 ' ^ ^ 



where the right hand side is clearly decreasing in n, I and N . 
Proof: Note that for n > 0, £ > and A^ > 3, 

V9A^2 _ + 65 + 4^2 + 4^ 
3A^ + 4n + 2£ - 5 

since then 



< 1 (8.14) 



(3A^ + 4n + 2^ - 5)2 - {m"^ - ASN + 65 + + Al) 

= {2AN - 40)n + {12N - 24)£ + 16n^ + 16n^ + 18A^ - 40 > . (8.15) 



□ 



8.4 LEMMA. For N > 4, g2iN) is a decreasing function of N . 
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Proof: Let h{x) = (1 - ^^'Z^, so that g2iN) = h{N — 1). Computing the derivative of 

ln(/i(x)), one finds that it is negative for x > 3. □ 

8.5 LEMMA. For n >0 and N >3, §3(11, N) is a decreasing function of n and N. 

Proof: For n a nonzero integer 

r(n + l)r(|iV-3) n n-1 1 



^.16) 



r(n + |7V-3) n + |A^-4 n + |iV-5 fiV - 3 ' 
Since each factor is less than 1 for > 3 and is a decreasing function of N the assertion follows. □ 

8.6 LEMMA. For N > 3, 94 (n + i, N) is a decreasing function of n + £ with 

lim g4{n + i,N) = 0. (8.17) 

n+£— >oo 

Moreover, forn + £>t = 3{N - 3) /2, 

(8.18) 

where 

/(^) = ^ — 4^ — ■ (^-'^^ 

Finally, for N > 5, {N — l)^f{N) is a decreasing function of N. 
Proof: Since 



r{n + £+^) r (n + £ + - I) _ (n + £ + |) 



r(n + ^ + |A^-|) r(n + £+|) n + £ + ^N ^ 

for > 3, it follows that for a fixed nonnegative integers A^, 

r{n + £ + l) 



5 < 1 (8-20) 



r{n + £ + lN-l] 



i.21) 



is a decreasing function of n + £. Hence, for n + £ > 3{N — 3) /2, 

T(n + £ + ^) _rf|iV-3) 



< 



r(n + £ + fA^-|) - r(3iV-7) 

This together with the definition of 5(4 proves the first inequality in (IS.lSp . Use of the duplication 
formula for the T function yields 

r(|iV-3)'_ V^r(|iV-3) _ ^(|Ar-4)r(|Ar-4) ^ V^(|iV-4) 



r(3iV-7) 23^-8r(|Af-|) 23JV-8r (lAT _ 7) - 23^-8 ■ 

This implies the second inequality in (jS.lSp . A check of the logarithmic derivative of {N — l)^f{N) 
shows it is negative for N < 5. 
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□ 

Now, combining the results in the last four lemmas, we have that for > 3 and n + ^ > = 
3(iV-3)/2, 

^ ^lA^) ^ '^'(^) (8-22) 

where 

= V ( 3N + 4n + 2e-5 + ^ 92iN)gsin, N)g,in,l, N) , (8.23) 



and 



.\N) = ^-^l^^^^ + l]g,{N)f{N), (8.24) 



where §2, gs and / are given by (|8.12|) and (|8.19|) . 

We are now ready to prove the main theorem of this section: 

Proof of Theorem 18.11 First, we take care of large values of N. By Lemmas 18.41 and 
(A^ — 1)^«;(A^) is a decreasing functions of N for N > 5. Direct computation shows that at = 12, 
this quantity is less than one. Hence for > 12, k,{N) < {N — 1)~^ = Kq2- For £ > £*, so that 
(j8.22p is satisfied, this proves (j8.ip for A^ > 12. On the other hand, if 2 < ^ < we have this from 
Lemma 18.21 or (|8.7p . Thus, in any case, (j8.ip is valid for A^ > 12. 

For 4 < A^ < 11, we again use Lemma 18.21 or (j8.7p for 2 < i < i* , and computation of kn/- 
By (|8.17p . for each such A^ there is a finite value A;(A^) so that we need only consider values of 
n + £ < k{N). Checking these cases, we obtain (j8.ip and (j8.2p . 

We finally turn to A^ = 3, which requires the greatest amount of computation. First for n = 0, 
we have from (18. 7p that 

Ko/(3) = 



^ 2 

so Ko,i(3) = -1/2 and |ko/(3)| < 1/3 for ^ > 2. 

The exact forms of the eigenvalues are simple enough to be useful for n = 1 and 2 as well. We 
have: 

. (m (-1)^+^ [2^Ar + 3(Ar - 1)] 
^uiN) - (AT _ iY+2 

and 

(-1)^((4£2 + IQi + 15)A^3 _ (8^2 ^ ^ QQ^j^2 ^ _^ ^Q^-^jy _ 12) 



K2/(iV) 



3(3Ar-4)(A^- 1)^+4 
Specializing to A^ = 3, 

.1,(3) = (-1)^+^^ 
so that 1^1, (3)1 < 3/8 for i >2. Likewise, for A^ = 3, 

£ (7 + 3£) , 
'^2/(3) = ^^^(-1) 

which implies that 

1^2/(3)1 < |k22(3)| = ^ (8.25) 
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For higher values of n, we estimate ^ by means of ^. Since i* = for = 3, we may use 
Lemma 18.61 for ah i, and then by (j8.17p . for each fixed n, there is a maximal value i{n) that need 
to be considered, and a maximum value of n that need to be considered. The following table gives 
the values of n,£ and k^£(3) when k^^(3) < 1/4. The monotonicity of k^£(3) in n and I shows 
that K^^(3) < «{3) for n > uq and i > io where (no,^o) is chosen from the table. □ 



n 


£ 




n 


£ 




n 


£ 




n 


£ 




3 


1253 


0.10562 


20 


210 


0.10547 


37 


90 


0.10543 


54 


34 


0.10514 


4 


989 


0.10562 


21 


198 


0.10559 


38 


86 


0.10531 


55 


31 


0.10538 


5 


817 


0.10556 


22 


188 


0.10546 


39 


82 


0.10527 


56 


29 


0.10506 


6 


694 


0.10561 


23 


178 


0.10552 


40 


78 


0.10528 


57 


26 


0.10538 


7 


604 


0.10555 


24 


169 


0.10549 


41 


74 


0.10537 


58 


24 


0.10511 


8 


533 


0.10560 


25 


161 


0.10537 


42 


70 


0.10551 


59 


21 


0.10551 


9 


477 


0.10558 


26 


153 


0.10540 


43 


67 


0.10523 


60 


19 


0.10529 


10 


431 


0.10558 


27 


145 


0.10558 


44 


63 


0.10552 


61 


17 


0.10509 


11 


393 


0.10555 


28 


138 


0.10561 


45 


60 


0.10534 


62 


14 


0.10560 


12 


360 


0.10561 


29 


132 


0.10542 


46 


57 


0.10521 


63 


12 


0.10545 


13 


333 


0.10548 


30 


126 


0.10534 


47 


54 


0.10512 


64 


10 


0.10534 


14 


308 


0.10558 


31 


120 


0.10540 


48 


50 


0.10562 


65 


8 


0.10523 


15 


287 


0.10554 


32 


114 


0.10554 


49 


48 


0.10508 


66 


6 


0.10514 


16 


268 


0.10556 


33 


109 


0.10543 


50 


45 


0.10512 


67 


4 


0.10509 


17 


251 


0.10558 


34 


104 


0.10540 


51 


42 


0.10521 


68 


2 


0.10506 


18 


236 


0.10554 


35 


99 


0.10546 


52 


39 


0.10535 


69 





0.10503 


19 


222 


0.10560 


36 


94 


0.10562 


53 


36 


0.10552 









The remaining values can be computed from the exact formula for Kn,£(3) from ()7.ip . and the 
results are all consistent with (j8.3p . □ 

9 The determination of the spectrum of P 

For given values of N ^ n and £, let iJLn,i{N) be the eigenvalue of P corresponding to the eigenvalue 
i^n/{N) of K through Theorem 18. 11 where we use (|5.7p if Kn/{N) > 0, and use (15. 8p if Kn/{N) < 0. 
(This is the relevant choice, as we are concerned with the largest eigenvalues of P.) 
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Consulting the calculations in fTTO]) for n + ^ = 1, and in frT2]) . ([7l3|) and (fTlij) forn + ^ = 2, 
and finally the bounds in Theorem lS.ll for n+i > 2, we see that for all > 3, the largest eigenvalues 
of K is and the least (most negative) is ko,i = '^i.o- Thus, turning to Lemma l5.ll and using 
the positve eigenvalue in (j5.7p . and the negative one in (jS.Sp . we see that the positive one yields 
the greater value for each N. Thus, the gap eigenvalue of P, fij\[, is given by 

, , 3iV - 1 , , 

/Xjv = m,i(iV) = 3(jy_-^^2 • (9-1) 

Use of this result in p.Sp would yield a strictly positive lower bound on A^r, uniform in N, but, 
as we have said above, it would not yield the exact lower bound. To obtain this, we now carry out 
the strategy outlined in the introduction. 

First, we combine Lemma l5.ll and Theorem 18.11 to produce the information necessary for the 
application of Lemma [4.1l We must now make a choice of the thresholds fi^ that appear in Lemma 
14. li The choice we shall make is based on trial function computations with Q that suggest that 
the gap eigenfunctions are the ones specified in Theorem I l.li 

N 

Notice that in Theorem ll.il the formula given for Am is of the form C for some constant 

' ^ N — 1 

C. This value can be guessed by computing the eigenvalues of Q on the invariant subspace of 
polynomials of degree 4 or less in the Vj . If we are to prove this guess correct using ()4.2p of Lemma 
14. H we require a value of /x^ such that 

N , ^ ,N -1 N , , 

[l-f^*N)T^>T^ , (9.2) 



iV-r "^"^ N -2 - N -1 

at least for > 4. (The guess is valid only for — 1 > 3. For — 1 = 2, there is a different value 
of A2 which has been determined already in Section 2.) 
The largest value of /i^ that will satisfy (19. 2p is 

This turns out to be an eigenvalue of P: Indeed, we have found in ()7.1ip that k;o,2 = 1/(A^ ~ 1)^- 
Furthermore, we have found in (jT.lOp that «:o,i = ki,o = — 1/(A^ — !)• Using first of these results in 
(|5.7p of Lemma 15. H and the second in (|5.8p we find 



^J'0,2 — Ati,o — A*o,i 



A^- 1 



For A^ = 3 we need to make a different choice, as the spectrum of Q is quite different for N = 2 
and for A^ > 3. The choice that shall work is fi^ = /i2,2(3). Since ^2,2(3) = 13/40, we have from 

LemmaOthat //2,2(3) = ^(1 + 2(13/40)) = ^. Thus, 

Now, to apply Lemma l4.ll we need the eigenspaces of P for the eigenvalues /i satisfying 1 > 
/i > n^. By Theorem 18.11 and (j7.13p . for A'' = 3 and A^ = 4, there is just one such eigenvalue, 
namely //i^i(4), the gap eigenvalue, and for N > 5, there are two: fii^i{N) and ^2,o(A^)- 
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Let En^£ be the eigenspace of P corresponding to the eigenvalue Hn,ei^)- For all values of 
n and i with n + £ < 2, we have determined the corresponding eigenfunctions of K, and thus, 
through Lemma 15.11 the corresponding eigenfunctions of P. Thus, we have the following explicit 
descriptions of the En/ for n + £ < 2: 

First, for n + £ = 1, the eigenvalues of K are negative, and so by Lemma [5 .11 the eigenfunctions 
are antisymmetric. If we are only concerned with the spectrum of Q on the subspace of symmetric 
functions (which is all that is of significance for Kac's application to the Boltzmann equation), 
we can ignore these eigenspaces. However, they turn out to be very simple. The n = 0,£ = 1 
eigenfuctions of K are degree one spherical harmonics, and the n = 1,1 = eigenfuctions of K are 
degree one Jacobi polynomials in Hence 

Eq^i is spanned by the functions vf — , a = 1, 2, 3 and i < j , (9-5) 

while 

Ei^Q is spanned by the functions — \vj\'^ ,^ < j , (9-6) 

Next, for n-\-£ = 2, the eigenvalues of K are positive, and so by Lemma lS.H the eigenfunctions 
are symmetric. The n = 0,1 = 2 eigenfuctions of K are degree two spherical harmonics, and so 
have the form 

3 

f0,2(.v) = Yl ^o.,pV^^^ 
o,/3=l 

for some traceless symmetric 3x3 matrix A. Hence, by Lemma |5.H 

N 

i?o,2 is spanned by the functions ^^/o,2('yj), i^-"^) 

with /o,2 given as above. 

For n = 1, £ = 1, the eigenfunctions of K are the product of a degree one spherical harmonic, 
and a degree one Jacobi polynomial in \v\'^ . When we sum over the particles, the constant term in 
the Jacobi polynomial drops out due to the momentum constraint, and we see that 

N 

Ei^i is spanned by the functions (9.8) 

where 

/i,i(^) = b|V a = 1,2,3. 

Finally, for n = 2, ^ = 0, the eigenfunction of is a degree two Jacobi polynomial in 
After summing on the particles, the linear term can be absorbed into the constant by the energy 
constraint, and so we see that 

N 

£^2,0 is spanned by the function ^^/2,o('yj)) (9.9) 

i=i 

where 

/2,o('y) = l^^l'* - / \v\'^duN . 

Jb 
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We close this section with a lemma that we shall need to prove Theorem 11.21 There we shall 
need to know the next largest eigenvalue of P below max„+^<2 f^n/i^)- One might guess that this 
occurs for some values of n and I with n + i = 3, but this is not the case: By (jS.Sp of Theorem l8.lt 
and (j7.12p . we se that for = 3, the most negative eigenvalue of K is —1/2, and by Lemma |5. 11 
this corresponds to the eigenvalue 1/2 of P. On the other hand, the largest eigenvalue of K appart 
from Ki^i(3) is ^2,2(3) = 13/40. This corresponds to the eigenvalue 11/20 of P. Since 11/20 > 1/2, 
we do indeed have that 

sup ^„,£(3) = /X2,2(3) = — . 

n+£>2 

It seems likely, on the basis of computations that we have made, that in fact 

sup ^Xn,iiN) = fX2,2{N) (9-10) 
n+e>2 

for all > 3. However, for the proof of Theorem 11.21 that we require is: 
9.1 LEMMA. For iV = 3,4,5,6 and 7, [9AU\) is true. 

Proof: Note that the case = 3 has already been proved in the remarks above. To deal with 
the other cases, we proceed essentially as in the proof of Theorem 18. H using ()8.23p to reduce the 
numbe of cases to be checked to a finite number, and then checking these. We will therefore be 
brief in our remarks on the remaining cases. 

Perhaps the most important point to recall is that (18.230 is valid for n + l> 3{N — 3)/2. since 
the right hand side evaluates to zero for = 3, we could use it without restriction. For N = 7 
though, 3(A^ — 3)/2 evaluates to 6, and so we may only use (|8.23p for n + £ > 6. So these cases 
must be checked be direct computation of the eigenvalues using [7?T|) . and then converting these to 
eigenvalues of P using Lemma |5. 11 

Then, using (|8.23p for n + £ > 6, one finds that 

f^lA'^) < '^1,2(7) 

unless < n < 6 and < £ < 27. Computing the rest of the eigenvalues of P in this 6 by 27 
rectangle, we find that stated result is true for N = 7. 

A similar analysis takes care of A = 4, A = 5, and A = 6. □ 

We shall not need to know the corresponding eigenfunctions in our application of Lemma 19. H 
since we will only be concerned with the eigenspaces of eigenvalues lying strictly above /i2,2(A), 
and those have been determined already in this section. 



10 The spectrum of Q on invariant subspaces containing 
eigenspaces of P 

For each n and £, let Vn/ be the smallest invariant subspace of Q constaining E^^t- As we shall 
see, for n + £ < 2, Vn/ = En/ except for n = 2, £ = 0, in which case V2,o is two dimensional, while 
Ei^Q is one dimensional. This is established in the next lemma, which also specifies the spectrum of 
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Ags = 1 - (1 - ^2) ,,,, , (10.3) 



Q on these invariant subspaces. The eigenvalues of course depend on the particular choice of b in 
the definition of Q, but in a very simple way: The dependence on b is only through the quantities 
(1 — Bi) and (1 — B2), where Bj is the jth moment of b, as defined in (jl.Sp . 

10.1 LEMMA. Every nonzero function in Eq^i and in Eifi is an eigenfunction of Q with eigen- 
value 

Ao^i = A?o = l-(l-^i)^ > (10.1) 

so that Vo,i = ^0,1 CLiT'd Vi^ = -Ei,o- 

Every nonzero function in Ei^i is an eigenfunction of Q with the eigenvalue 

A?i = l-(l-i?2)^^^^ , (10.2) 

so that Vi^i = Ei^i. 

Furthermore, every nonzero function in £'0,2 is an eigenfunction of Q with the eigenvalue 

3 

'2(iV"- 1) 
so that Vb,2 = -E'o,2- 

Finally, while 1^2,0 is larger than £^2,0 > there are only two eigenvalues of Q in V2,o- These are 

l-{l-B2) ( ^^^^_^^ [(2iV - 1) ± ViV2_3^r + i] ^ . (10.4) 
For all N > 3, the largest of these eigenvalues is Xfi . 

Before beginning the proof, we note that if (l/2)6(x)dx is a Dirac mass at x = 1, the collisions 
are all trivial (zero scattering angle), and thus Q = / in this case. But also in this case {1 — Bi) = 
(1 — B2) = 0, so all of the eigenvalues A^^ listed above are 1 — as they must be for Q = I. 

Proof: We begin with the last case, case n = 2, i = 0, which is the most involved. Consider the 
function 

N 

4> = ^\vi\' (10.5) 

i=l 

and note that cf) — J"^^ (pdajy spans £2,0) ^ we have noted above. 

One simple way to calculate Q(p is to take advantage of the permutation symmetry of Q: Define 
the symmetrization operator S by 



Sf{vi,. . . ,vn) = -^^/(u^(i),...,?;^( 



where the sum runs over all permutations vr of {1, . . . , N}. Then it is easy to see that 

5(hh = ^<A. 

Thus, since SQ = QS, 

Q<P = NSQ{\vi\^) . 

One now directly calculates Qdt'il'') and then symmetrizes. In carrying out the calculation, we 
make use of: 
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10.2 LEMMA. Let c and d be any two vectors in "M? , and let e be any unit vector in . Then 
with Bi and B2 defined as in il.5\) . we have the following identities: 



/ (c • a)b{e ■ a)da = (c • e)Bi 



and 



-j^ 35 1 

(c ■ a){d ■ a)b{e ■ cr)dcj = (c • d) — h (e • c)(e • d) 

52 2 2 



Proof: We choose coordinates in which c and e span the x, z plane with e 



" " 




" ci ■ 





and c = 





1 




c3 



, the computations are easily accomplished. 



□. 



sm sm (p 
Then with a = cos 6 sin 99 
cos 6 

Now to compute Q(j), go back to the definition of Q given in (jl.4p : and note first of all that with 
Tjiv) = \vi\^, unless i = 1, 

v{Ri,j,a{v)) = r]{v) . 

Hence 



da . 



Then from (II. 
r]{Rij,aiv)) = 



Vl+Vj \vi-Vj\ 

H —a 



= ^ IbiP + biP + 1^1 - Vj\{vi + Vj) ■ af 
o 

= I ((bil^ + bjf )^ + 2(biP + hjf - Vjlivi + Vj) -a + lvi- VjlHivi + Vj) ■ af) . 

(10.6) 

Integrating over S"^ using Lemma 110.21 vields 



/ rj{R,j,Av))b 



\Vi - Vj\ 



2\2 



:i\vi\' + \v,\') 



+ -(\v,\^ + \v/)B,{\v,\'-\v,\') 

I 1 2 1 1 2 ^ — / / \ / \ \ 2 — 1 

+ \Vl-Vj\ \Vl+Vj\ — \- {{Vl - Vj) ■ [Vl +Vj)) 



16 
(10.7) 



The right hand side simplifies to 

^(\v,\' + \v,\^ + 2\v^\'\v,\^) + ^{\v^\'-\v,\^) 



+ {\vi\^ + \v,\^ + 2\vi\''\v,\''-iiv,.v,f)^ 



16 



+ {\v,\' + \v,\'-2\v^\'\v,\'f-^ 



(10.8) 
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It is now a simple matter to carry out the sum on j > 2. Using the identities 



N 

^ \Vj\^ = (t){v) - \vi\ 

i=2 

and the symmetrizing, one finds 



I-B2 

N 



and 



N 

El 

J=2 



(v) - \vi\ 



+ — ttV' 



N-r (iV-l) (N-l) 

where 

For > 4, the two functions </) and ^ are linearly independent, although for = 3 they are 
not. In fact for = 3, one has the identity 



(10.9) 
(10.10) 



flO.lll 



Evidently, for > 4, 1?!) — (^do"Ar is not an eigenfunction of Q, so that £^0,2 is not an eigenspace 
of Q. We are required to compute QV- 

We again take advantage of the permutation symmetry, and note that 



S{{v\ ■ V2) 



^ V and Q^=^^SE_^S((^{y^-V2f) 



N{N - 1) " " 2 

We carry out the calculation in the same way that we calculated Q(p, and find 



1-B2 

N 



3V^ + 



(N-l)' 



1 

N 



fl0.12) 



We see that the subspace spanned hy (j) — f-^^ (j)daN and ip — J^^ (j)da]\f is invariant under Q. 
Using (jl0.12p and (jlO.Op we easily find that the two eigenvalues of N{I — Q) on the two dimensional 
space V2,o are the eigenvalues of 



which are 



I-B2 
N-l 

I-B2 



N + 1 1 
Af - 3 3A^ - 3 



{2N - 1) ± \/Ar2 - 3A^ + 1 



N-l 

The minus sign clearly gives the lesser of these, and gives the gap for N{I — Q) on V2,o- From here, 
one easily deduces ()10.4p . 

A further, much simpler calculation shows that the three functions 



TV 



^1,1 = ^\^k 



2„,a 



k=l 



where a indexes the components, are also eigenfunctions of Q, more precisely 

I-B2 



1 



A^- 1 



(10.13) 



(10.14) 
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Thus the unique eigenvalue of A on Fi^i is 

Q _ I-B2 
\i - 1 aTTT • 

For -Eo,2, a simple computation shows that the functions 

k 

where a ^ (3 are indices for the components, are also eigenfunctions for Q, in fact 

2{N - 1) 



Thus, Vb,2 = -^0,2; and the unique eigenvalue of Q on this subspace is 

- 1 - ^ 

An o — J- — 



Finally, we consider the spectrum on Q on the eigenspaces of P corresponding to n + ^ = 1. In 
this CctSG, clS noted above, the eigenfunctions are antisymmetric, so that if we are only concerned 
with the spectrum of Q on the subspace of symmetric functions (which is all that is of significance 
for Kac's application to the Boltzmann equation), we can ignore these eigenspaces. However, if we 
define ??o,i(^^) = vi — V2 and riifi{v) = \vi\'^ — |i'2p, we find, as above, that 

QVifi = (1 - (1 - Bi)j^^ Vi,o and Q%,i = (l - (1 - %,i . (10.17) 

□ 

Now that we have all of our eigenvalues we need to order them. By a simple comparison, we 
determine that for all N, the largest eigenvalue of Q on our three invariant subspaccs with n + £ = 2 
is Xii- This is true for all choices of 6, since the only dependence on b in these eigenvalues is a 
common factor of (1 — B2). 

It is worth noting, however, that for large N, 

A?, = l-(l-S2)l + o(^) and A?o = l-(l-i?2)^ + o(^) , 
so that these eigenvalues merge as N tends to infinity. Still, for all finite N, 

^?i = l-(l-^^)]^ 



is strictly larger. 

Next, the invariant subspaces of Q with n + £ = 1 are also eigenspaces of Q with the eigenvalue 

A?i = A?o = l-(l-5i)^- 



In summary, the largest eigenvalue of Q on the invariant subspaces Vn/ in L^(X7v, dcrjv) with 
n + £ = l,2 and N > 3 is either 

l-(l-52)— ^ or 1-{1-Bi ^ 



N -1 ^ ' 



CGL May, 2007 



36 



depending on which of these is larger. In particular, 

A3 <min{(l-B2),(l-5i)}^ (10.18) 

With the above arguments we have all the ingredients needed to prove Theorem ll.lt 

Proof of Theorem ll.lt First, we wish to apply Lemma |4. II to estimate A3 in terms of A2. In 
(j9.4p . we have set fi^ = 11/20, and with this choice of the threshold, we have seen that there is just 
one eigenvalue of P between /Xg and 1, namely the gap eigenvalue fi^ = /ii^i(3) = /Uo,i(3) = fiifi{3). 
Thus, from Lemma l4.1l and the eigenvalue computations in Lemma llO. It either the gap eigenvalue 
of Q for iV = 3 is 

max I l-(l-i?2)^^ , l-{l-Bi)j^ I , (10.19) 

or else 

A3 > - f 1 - — ^ A2 = — A2 . (10.20) 
2 V 20y 40 ^ ^ 

If (|10.19|) does give the gap eigenvalue of Q for = 3, then 

A3 = min{ (1 - Si) , (1 - S2) }^ . (10.21) 

Since according to Lemma |4.H at least one of (|10.20|) and (|10.21|) is true, the condition 

27 3 

-A2 > - min{ (1 - B,) , {I-B2)}, (10.22) 

and (110. ISp ensure that (110. 2ip is true, and thus gives us the gap eigenvalue for = 3. Note that 
the condition (jl0.22p is equivalent to the condition (jl.Gp in Theorem 11.11 
Now we proceed by induction. For any n > 4, assume that 

N - I 

An-1 = min{ (1 - Bi) , (1 - B2) . (10.23) 

In (j9.3p we have set 

^ _ 1 

for all A'^ > 4, and we have seen that the only eigenvalues ^ oi P with 1 > /U > (jl*^ are the gap 
eigenvalue = ^i^i{N) = ij,q^i{N) = fiifi{N), and for > 4, fi2fl{N). Thus, by Lemma ITT] and 
the eigenvalue computations in Lemma llO.H either the gap eigenvalue of Q for A" is 

max I 1 - (1 - B2)j;^ , 1 - (1 - Bi)j^ | , (10.24) 

or else 

Aa. > ^ (1 - ^) A^_i . (10.25) 

There is strict inequality in ()10.25p since all remaining eigenvalues of P not taken into account in 
(|10.24p are strictly less than /u^. By the inductive hypothesis (|10.23p . (|10.23p yields 

Ajv>min{ (l-Bi), {l-B2)}j^. 
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This is impossible, as the trial functions leading to (jl0.24p yield the upper bound 

N 

A7V<min{ (1-Si) , (1-^2) }^^-^ . (10.26) 

Thus equality holds in (110. 26p . which completes the proof of the inductive step. Because of the strict 
inequality in ()10.25p . the only eigenfunctions with the gap eigenvalue are found in the invariant 
subspaces considered here; i.e., in the Vn/ with < n + i < 2. By the results of Lemma ri0.1|, this 
yields the statement in Theorem 11.11 concerning the gap eigenfunctions of Q. □ 

Proof of Theorem II. 2t We proceed as in the previous proof, except that for low values of N, we 
shall use a different choice for the threhold /x^, namely 

f^N = f^2,2{N) . (10.27) 

We know from Lemma |9. II that for all N < 7, 

A*n,K^) < /^2,2(A^) for ah n + £>2 . 

Thus, at least for such N, the only eigenvalues fi oi P with /i > /i^ = fj,2,2i^) are those with 
n + i < 2. We have already computed the gap for Q on the invariant subspaces containing these 
eigenvalues, and we have found that the gap in these subspaces is 

N 



min{ (1-Si) , (1-^2) } 



iV - 1 



If for any Nq > 3, it turns out A^Vq = Atvq, the gap on the whole space, then we can switch 
from that point onwards to the use of //^ = Ho,2{^) as in the proof of Theorem 11.11 to show that 
Atv = Aat for all > A'^O) and that the eigenfunctions are exactly as claimed for all > A^o- 

We now show that it is always the case that Ajvq = ^Nq for some A^'o < 7. To do this, pick any 
value A''i > 4, and suppose that for 3 < j < A?"!, we have 

A,- <min{ (l-5i) , (I-S2) • (10.28) 

J - 1 

Then by Lemma l4.1l and Lemma l9.ll using the value fj,^ = H2,2{j), for 3 < j < A^^i, we have 

N 

AjV, > ^n(l-/^2,2(i))A2 . 

J=3 

Using the hypothesis A2 = 2(1 — Bi), we have 



>-in(i-M2,2(j))2(i-i?i 



2 

i=3 



Of course we can rewrite this as 



^ f (&i-m2(.o) 1(1 -.2,2(3)) (n 



^1-1 1/^4(1-/^0,2(7) ) 10 



-^(l-i?l) , 

(10.29) 
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since, as in the last proof, (1 — /X2,2(3)) = 9/20, and 

^11(1 -/^0,2U)) 



J =4 



Now, by direct computation, we find 



TT- (1 - /i2,2(j) _ 558018643 10 
fJl (1 - /^o,2(j) ~ 495720000 ^ ¥ " 

For A''i > 7, this would lead to 

Aiv.>^(l-5i), 

and this is impossible, since we have a trial function showing that the gap cannot be so large. 
Hence it must be that (110.28P is false for some j < 7. By what we have said above, from this point 
onward, we can proceed as in the proof of Theorem ll.il and we obtain Theorem 11.21 □ 

While the results presented here cover a very wide range of models, it is possible to come up 
with choices of b for which A2 7^ 2(1 — Bi). If one found a need to deal with such an example, 
one might have to go deeper into the spectrum of P. It is very likely that Lemma 19.11 holds for 
all > 3, based on extensive computation. These computations also show that as A'^ increases, 
/i2,i(-/V) comes very close to fJ-2,2{N), so that to get much more leverage, one would need to compute 
all of the eigenvalues of Q on the smallest invariant subspaces of Q that contains both of these 
eigenspaces of P. This could be done using the methods presented here, but the computations 
would be considerably more involved than the ones we have presented in this section. Thus, having 
treated a wide range of models, we shall conclude our discussion of Q here. In the brief final section, 
we discuss a point we raised earlier concerning bounds on Jacobi polynomials. 



11 Bounds on Jacobi polynomials 



As alluded to in Section 8 the identity (j7.9p . together with the trivial bound on the |Kn/| < li 
which comes from the fact that if is a Markov operator, will for certain regions provide a stronger 
bound than (jS.Sp . the bound of Nevai, Erdelyi and Magnus. We close this section by showing how 
(|7.9p can be used to obtain better bounds. 
To begin, write 

^2/3-1 ^n'^(-l + 2fe') < 2e r(n + l) 2 + ^^^TW r(n + /3 + l) r(a + 1)^ 

Pn'^il) ~ vr 6(l-62)a+i/2 2n + a + /3 + ir(ra + a + /3 + l)r(n + a + l)' ^ 

where (3 = 1 + 1/2 with I an integer. In regions where the right hand side of the above equation 
becomes larger than one the simple bound 

Pn'^il) 
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becomes stronger. In the region 2n + 1 < a < /3, we find 2n+Q+/3+ T ^ i- This plus Stirling's 
formula with the remainder yields, 



2e r(n + l) 2 + vV + /32 r(n + /3 + l) r(a + 1)^ 
V 6(1 - 62)a+i/2 2n + a + /? + 1 r(n + a + /? + 1) r(n + a + 1) 

e" r(n + l) a"+i/2-"/?-° (1 + 

^ V2^6(l - 62)^+1/2 (1 + ^ H±l )n+a+l/2 (l + B±l )/3 * 

^ e" r(n + 1) a°+i/2-"/3-« ^ 



2:^5(1 _ 52)0+1/2 22n+2a+2/3+3/2 ' 

where r = ( 1 — — — — — | f 1 — — — — I and n is assumed to be fixed. Choosing 

V 12(n + Q + /3 + l)y V 12(n + a + l)y 

6(1 — 6^)"+^/^ so that the last inequality is greater than one provides a region where (j7.9p and 
l/^n/l ^ 1 does better than (j7.9p . It would be interesting to obtain better bounds on \Kn/\ by direct 
analysis of K, and to use these to sharpen the argument just made. 
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